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Abstract 


A  manufactured  product  often  possesses  residual  texture  which  was  either  inciden¬ 
tally  or  deliberately  acquired  during  its  processing  history.  This  is  particularly  true  for 
the  case  of  polymers,  where  the  ability  to  easily  preferentially  pre-orient  the  material 
in  specific  directions  is  exploited  in  order  to  obtain  a  higher  strength  product.  Specific 
examples  include  synthetic  fibers,  and  biaxial ly-oriented  films  and  containers.  The  re¬ 
sponse  of  the  pre-oriented/textured  product  to  normal  service  life  loading  conditions 
will  differ  considerably  from  that  of  a  product  composed  of  isotropic  material.  This 
paper  addresses  the  issue  of  the  effects  of  texture  on  the  deformation  behavior  of  glassy 
polymers.  Here,  the  physically-based  constitutive  model  of  Boyce,  Parks,  and  Argon 
[1987]  describing  the  rate,  temperature,  and  pressure  dependent  inelastic  deformation 
of  initially  isotropic  glassy  polymers  is  extended  to  include  effects  of  pre-orientation, 
i.e.  initial  texture,  via  the  use  of  appropriate  initial  conditions  on  internal  state  vari¬ 
ables.  The  modified  model  is  then  utilized  in  an  analysis  of  the  effects  of  texture  on 
the  yield  of  glassy  polymers  and  the  shear  localization  which  normally  follows  yielding 
in  oriented  polymers.  These  results  are  compared  with  trends  found  in  experiments  as 
reported  in  the  literature.  The  effectiveness  of  the  proposed  model  for  the  present  appli¬ 
cation  is  also  compared  with  earlier  models  of  yielding  of  anisotropic  materials  such  as 


Hill’s  criterion. 
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19.  Abstract  (Continuation) 


The  modified  model  is  then  utilized  in  an  analysis  of  the  effects  of  texture  on 
the  yield  of  glassy  polymers  and  the  shear  localization  which  normally  follows 
yielding  in  oriented  polymers.  These  results  are  compared  with  trends  found  in  ex¬ 
periments,  as  reported  in  the  literature.  The  effectiveness  of  the  proposed  model 
for  the  present  application  is  also  compared  with  earlier  models  of  yielding  of 
anisotropic  materials  such  as  Hill's  criterion. 
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1  Introduction 

Texture  affects  both  the  elastic  and  the  inelastic  behavior  of  a  polymer.  However,  the 
effect  of  texture  on  the  inelastic  response  is  far  greater  than  that  on  the  elastic  response, 
especially  for  glassy  polymers.  Hadley  [1975]  has  summarized  the  effect  of  orientation 
on  the  small  strain  elastic  properties  of  polymers.  He  reports  that,  whereas  the  elastic 
tensile  modulus  of  some  semi-crystalline  polymers  may  increase  by  more  than  an  order 
of  magnitude  at  large  draw  ratios,  the  largest  increase  in  the  modulus  of  amorphous 
polymers  was  found  in  PMMA  and  was  only  30%  (Wright,  et  of.  [1971]).  This  increase 
in  elastic  modulus  with  orientation  is  very  small  and  would  have  little  impact  on  the 
material  undergoing  large  inelastic  deformation,  where  the  effect  of  orientation  on  the 
yield  strength  of  the  material  is  considerable.  Therefore,  for  the  purposes  of  this  work, 
orientation  effects  on  the  elastic  response  of  the  material  will  be  neglected. 

The  effect  of  texture  on  the  inelastic  response  of  glassy  polymers  will  now  be  exam¬ 
ined.  In  order  to  best  understand  the  oriented  polymer,  the  primary  physical  mecha¬ 
nisms  governing  plastic  flow  in  initially  isotropic  glassy  polymers  will  first  be  reviewed. 
We  will  then  present  a  three-dimensional  constitutive  model  for  initially  isotropic  glassy 
polymers  proposed  by  Boyce,  Parks,  and  Argon  [1987a].  The  presence  of  an  initial 
texture  on  the  basic  macromolecular  structure  of  the  polymer  and  its  impact  on  the 
deformation  resistances  will  then  be  discussed.  These  observations  will  lead  us  to  in¬ 
corporate  texture  effects  into  the  constitutive  model  via  the  use  of  initial  conditions 
on  internal  state  variables.  Finally,  the  plastic  flow  in  oriented  glassy  polymers  will  be 
quantitatively  examined  using  the  model  in  conjunction  with  the  finite  element  method. 
In  particular,  the  effects  of  texture  on  the  initial  yield  and  subsequent  shear  localization 
will  be  analyzed.  The  numerical  results  will  be  compared  with  behavior  trends  which 
have  been  observed  experimentally  and  reported  in  the  polymer  literature. 
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2  Inelastic  Deformation  of  Glassy  Polymers 


In  general,  the  yield  of  initially  isotropic  glassy  polymers  has  been  found  to  depend 
on  strain  rate,  pressure  and  temperature.  After  yielding,  polymers  often  strain  soften 
and  subsequently  strain  harden.  These  characteristics  of  yielding  in  glassy  polymers  are 
considered  to  result  from  two  physically  distinct  sources  of  internal  resistance  (Haward, 
Thackray  [1968],  Argon  [1973],  Ward  [1984]).  To  best  describe  these  mechanisms,  we 
assume  an  idealized  material  structure  consisting  of  a  random  array  of  long  molecular 
chains  which  intersect  one  another,  forming  points  of  entanglement  [Figure  1],  This 
results  in  an  isotropic  network-like  structure  much  like  that  of  rubber,  but  with  the 
chemical  crosslinks  replaced  by  physical  entanglements.  This  structure  is  now  loaded. 
At  yield,  an  isotropic  intermolecular  resistance  to  chain  segment  rotation  must  be 
exceeded;  this  is  a  thermally-activated  process.  Once  the  material  is  free  to  flow, 
molecular  alignment  occurs,  resulting  in  an  anisotropic  internal  resistance  to  further 
inelastic  deformation.  Below,  we  discuss  these  deformation  resistances  to  chain  segment 
rotation  and  subsequent  molecular  alignment. 

2.1  Intermolecular  Resistance 

The  intermolecular  resistance  to  plastic  flow  is  due  to  the  impedence  imposed  by  neigh¬ 
boring  chains  on  the  ability  of  a  chain  segment  to  rotate  either  individually  or  in  a 
cluster.  Argon  [1973]  has  developed  an  expression  for  the  plastic  shear  strain  rate, 
which  results  once  the  free  energy  barrier  to  chain  segment  rotation  in  an  equivalent 
isotropic  elastic  medium  has  been  overcome: 

y  =  if0exp[-^(l-(l)l)];  (1) 

where  -yo  is  the  pre-exponential  factor;  A  is  proportional  to  the  activation  volume/ Boltz¬ 
mann’s  constant;  6  is  the  absolute  temperature;  sq  =  is  the  athermal  shear 

strength,  n  is  the  elastic  shear  modulus,  and  u  is  Poisson’s  ratio;  r  is  the  applied  shear 
stress.  This  expression  has  been  extended  by  Boyce,  et  al.  [1987a]  to  include  the  effects 
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of  pressure  and  strain  softening.  In  the  modification,  the  athermal  shear  strength  sQ  is 

t 

replaced  by  s  +  ap,  where  p  is  the  pressure,  a  is  the  pressure  dependence  coefficient,  \ 

and  3  is  now  taken  to  evolve  with  plastic  straining.  The  initial  condition  on  3  is  the 
annealed  value  s0.  The  shear  resistance  is  phenomenologically  modelled  to  decrease 
with  plastic  straining  until  reaching  a  “preferred”  structure,  represented  by  s„,  via 
3  =  h(l  —  where  h  is  the  rate  of  resistance  drop  with  respect  to  the  plastic 

strain,  and  sM  generally  depends  on  6  and  ip.  The  evolution  of  s  with  plastic  strain 
models  the  macroscopic  response  of  true  strain  softening.  Softening  is  considered  to  be 
the  result  of  local  dilatation  which  accompanies  micro-shear  banding  during  plastic  flow 
(Argon,  tt  al.  [1985]).  Any  fluctuations  in  local  volume  would  obviously  impact  the 
intermolecular  resistance.  This  leads  to  the  phenomenological  model  of  s  given  above. 

Thus,  3  is  taken  to  be  am  internal  state  variable  which  monitors  the  isotropic  resistance 
to  deformation.  The  shear  resistance  is  also  known  to  increase  with  pressure,  time,  and 
temperature,  i.e.  the  material  physically  ages  (Kramer  [1970],  Struik  [1978],  Haward 
[1980]).  Although  aging  is  not  currently  included  in  the  model,  it  is  important  to 
mention  this  phenomenon  prior  to  the  later  discussion  on  initial  conditions  for  internal 
state  variables. 

2.2  Network  Resistance 

Once  the  material  is  stressed  to  the  point  of  overcoming  intermolecular  barriers  to  chain 
motion,  the  molecular  chains  will  tend  to  align  along  the  direction  of  principal  plastic 
stretch.  This  action  decreases  the  configurational  entropy  of  the  system  which,  in  turn, 
creates  an  internal  network  stress  state.  Following  mechanics  terminology,  we  call  this 
internal  stress  the  back  stress  tensor,  B.  Due  to  the  rubbery  network-like  response 
of  glassy  polymers  to  plastic  deformation  (see  (Boyce,  et  al.  [1987a],  Ward  [1984])), 
the  back  stress  has  been  modelled  using  a  statistical  mechanics  of  rubber  elasticity 
model.  We  have  assumed  that  the  initial  network  structure  of  the  polymer  is  isotropic. 

Subsequent  network  orientation  is  taken  to  coincide  with  plastic  stretching,  and  B  is 


4 


A  »jt ^ ; *|V  *^‘«|>» A*|L*|tt*fcU%|,«il,>4i«*tl»> A^tC*^*^**^* W l*t' 


taken  to  be  coaxial  with  the  plastic  stretch.  In  terms  of  principal  components,  B  may 
be  expressed  as: 

B,  =  Bi{  Aj/Ax);  (2) 

where  D,  is  a  principal  component  of  the  back  stress  tensor;  Ay  are  the  principal  network 
(and  plastic)  stretch  components;  A i  is  the  locking  network  stretch  or  natural  draw  ratio 
of  the  material.  The  model  is  based  on  a  Langevin  spring  formulation  such  that  £, 
becomes  infinitely  large  as  Af/Aj,  approaches  1,  and  B<  =  0  for  the  isotropic  case  of 

a;  =  i. 

Above  we  gave  a  brief  summary  of  the  resistances  to  plastic  flow  in  an  essentially 
one-dimensional  framework.  These  ideas  have  been  incorporated  into  a  fully  three- 
dimensional  representation  of  inelastic  behavior  which  is  summarized  immediately  be- 


2.3  Three-Dimensional  Representation 

We  first  consider  the  deformation  of  an  initially  isotropic  body,  Bo.  The  body  is  loaded 
to  a  state,  Bt,  and  its  deformed  configuration  may  be  described  by  the  deformation 
gradient,  F,  at  this  time.  The  deformation  gradient  may  be  multiplicatively  decom¬ 
posed  into  elastic  and  plastic  components,  F  =  F*FP  (Lee  [1969]).  Here,  with  no  loss 
of  generality  (Boyce,  et  al.  [1987c]),  we  take  F*  to  be  symmetric,  i.e.  Fp  represents  the 
relaxed  configuration  obtained  by  elastically  unloading  without  rotation  (in  the  polar 
decomposition  sense)  to  a  stress-free  state.  The  rate  quantities  corresponding  to  this 
formulation  consist  of  the  velocity  gradient: 

L  =  FF’1  =  D  +  W  =  F*F*“l  +  F*FpFp-1F*~l,  (3) 

where  D  is  the  rate  of  deformation,  and  W  is  the  spin.  The  velocity  gradient  of  the 
relaxed  configuration  is  given  by  1/  =  FPFP_1  =  D*  +  Wp.  As  discussed  in  Onat 
[1987],  the  spin  of  the  relaxed  configuration,  Wp,  is  algebraically  defined  as  a  result  of 
the  symmetry  imposed  on  the  elastic  deformation  gradient.  The  rate  of  shape  change, 
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Dp,  must  be  constitutively  prescribed.  The  magnitude  of  Dp  is  given  by  the  plastic 
shear  strain  rate,  7P,  and  the  tensor  direction  of  Dp  is  specified  by  N  = 
normalized  deviatoric  portion  of  the  driving  stress  state,  T*,  at  a  material  point.  This 
driving  stress  state,  T*,  is  given  by  the  tensor: 

T*  =  T  -  ^F*BF*,  (4) 

«/ 

where:  T  is  the  Cauchy  stress  tensor;  J  is  the  volume  change  given  by  detF*;  and  B  is 
the  back  stress  tensor  due  to  strain  hardening  resulting  from  molecular  alignment.  We 
refer  to  T*  as  the  driving  stress  state  because  it  is  only  this  portion  of  the  stress  which 
continues  to  activate  plastic  flow.  The  back  stress  portion  is  an  internal  stress  which 
acts  as  a  resistance.  The  back  stress  tensor  is  taken  to  be  an  isotropic  function  of  the 
left  plastic  stretch  tensor: 

F p  —  VPRP;  polar  decomposition  (5) 

Vp  =  RprVpRp;  eigen  value  decomposition  (6) 

Bi  =  Bx(V’I\l).  (7) 

The  effective  equivalent  shear  stress,  r,  is  given  by: 

r  =  [IT“.T“]i.  (8) 

The  plastic  shear  strain  rate,  7P,  was  constitutively  prescribed  in  equation  (1).  However, 
r  is  now  taken  to  be  the  effective  equivalent  shear  stress  given  above  in  equation  (8). 

2.4  Effect  of  Initial  Texture  on  Micromechanisms  of  Flow 

Now,  we  consider  the  effects  of  pre-orientation  on  the  two  resistances  to  plastic  flow  - 
intermolecular  interaction  and  molecular  alignment.  The  intermolecular  resistance  is 
a  restriction  on  molecular  motion  due  to  surrounding  chains.  One  could  argue  that 
this  resistance  may  change  with  materia]  orientation  for  a  number  of  reasons.  At  very 
large  deformations,  after  much  molecular  alignment  has  occurred,  the  volume  fraction 
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of  sites  which  have  yet  to  overcome  intermolecular  barriers  has  greatly  decreased.  The 
presumably  higher  barriers  of  the  remaining  population  would  result  in  an  isotropic 
hardening  component  which  could  be  phenomenologically  accounted  for  within  the 
current  framework  of  the  model  by  an  increasing  athermal  shear  resistance  s  as  very 
large  deformations  are  reached.  We  also  note  that  the  intermolecular  resistance  is  a 
function  of  the  elastic  shear  modulus.  In  section  one,  we  concluded  that  the  effects 
of  pre-orientation  on  the  elastic  properties  are  small  and,  in  fact,  negligible  in  the 
current  context  of  application  of  this  work.  In  order  to  accurately  assess  the  effects 
of  developing  anisotropy  of  the  elastic  moduli  on  s,  an  equation  for  the  plastic  shear 
strain  rate  (equation  (1))  would  have  to  be  derived  considering  chain  segment  rotation 
in  an  anisotropic  elastic  medium.  This  effect  would  become  noticeable  only  at  very 
large  deformations,  at  which  point,  it  would  also  be  comparatively  small  with  respect 
to  the  asymptotically  increasing  back  stress.  Therefore,  in  our  modelling,  the  effects  of 
pre-orientation  on  the  intermolecular  resistance  will  be  neglected.  On  the  other  hand, 
the  network  resistance,  i.t.  the  back  stress  tensor,  is  obviously  greatly  affected  by  the 
presence  of  an  initial  orientation,  since  it  is  a  direct  result  of  molecular  alignment. 

The  above  discussion  leads  us  to  model  texture/orientation  in  glassy  polymers  as 
am  initial  network  stretch  in  the  material,  which  has  a  corresponding  initial  back  stress. 
In  other  words,  before  plastic  flow  commences  in  a  textured  polymer,  in  addition  to  the 
intermolecular  resistance,  an  initial  back  stress  must  also  be  overcome.  Figure  2  depicts 
a  schematic  of  the  true  stress-stretch  curve  of  an  initially  isotropic  glassy  polymer  which 
is  loaded,  unloaded,  and  then  reloaded  without  allowing  for  aging  to  occur.  The  effect  of 
the  initial  loading  (which  orients  the  material)  upon  the  stress  level  at  which  subsequent 
plastic  flow  commences  is  apparent  upon  reloading,  i.e.,  elevated  levels  of  “peak  yield” 
stresses  are  observed  in  the  direction  of  the  principal  orientation.  It  is  an  oriented  state 
such  as  this  that  is  now  our  starting  point.  Experiments  by  Ward  [1984]  on  cold- rolled 
and  spun  PET  and  by  Bahadur  [1973]  on  cold-rolled  PC  greatly  support  this  concept 
of  principal  orientation  as  a  description  of  state  in  glassy  polymers. 


2.4.1  Incorporation  of  Initial  Texture  into  Constitutive  Model 


We  now  consider  the  incorporation  of  an  initial  texture  into  the  above  constitutive 
model.  Prior  to  doing  so,  it  is  important  to  recall  the  idea  of  monitoring  and  describing 
the  “state”  of  any  material  by  the  use  of  internal  variables.  The  concept  of  internal 
state  variables  has  achieved  growing  acceptance  in  the  plasticity  field  in  recent  years 
(Kocks,  Argon,  and  Ashby  [1975],  Rice  [1975],  Onat  [1982],  Nemat-Nasser  [1983],  Loret 
[1983],  Dafalias  [1985],  Anand  [1985]).  In  principle,  an  adequate  number  of  internal 
variables  will  uniquely  describe  the  current  state  of  a  material  without  regard  to  prior 
history  of  loading.  Therefore,  if  given  the  instantaneous  values  of  the  relevant  internal 
state  variables  in  a  given  material,  one  should,  using  suitable  evolution  equations,  be 
able  to  accurately  predict  the  material  response  with  no  other  knowledge  of  the  process 
or  history  by  which  this  state  was  achieved. 

In  this  work,  we  are  putting  the  concept  of  the  internal  state  variable  to  full  use.  In 
the  model  for  the  initially  isotropic  glassy  polymer,  it  was  proposed  that  the  current 
state  of  the  material  may  be  described  by  the  stress  tensor,  T,  the  absolute  tempera¬ 
ture,  6,  the  athermal  shear  resistance,  s,  and  the  back  stress  tensor,  B.  Furthermore, 
the  back  stress  tensor  is  taken  to  be  a  state  function  of  the  plastic  stretch  measured 
with  respect  to  the  configuration  possessing  an  isotropic  network.  More  fundamentally, 
this  portion  of  state  description  is  described  by  a  network  stretch  tensor,  relative  to  an 
initial  isotropic  network  configuration.  The  coincidence  of  the  “plastic”  and  “network” 
stretches  follows  from  the  affine  connection  between  molecular  alignment  and  defor¬ 
mation.  For  the  case  of  the  oriented  glassy  polymer,  we  propose  that  the  constitutive 
model  for  the  initially  isotropic  glassy  polymer  also  properly  models  the  behavior  of 
the  corresponding  oriented  material  when  appropriate  initial  conditions  are  used  on 
the  athermal  shear  resistance,  s,  the  network  stretches,  which  will  be  called  A',  and,  of 
course,  if  any  residual  stresses  exist,  the  stress  tensor,  T.  In  turn,  the  initial  value  of 
the  back  stress  tensor  B  will  be  taken  as  a  function  of  (cf.  equations  5-7)  of  A*.  Below, 
we  will  not  consider  residual  stresses,  but  this  would  not  significantly  alter  the  basic 


frair  -  work  of  modifications  to  be  presented. 

In  order  to  best  describe  the  incorporation  of  preorientation  into  the  model  by 
the  use  of  initial  conditions  on  the  internal  variables  s  and  A',  we  will  consider  the 
deformation  of  two  bodies  —  the  initially  isotropic  body  and  the  initially  oriented  body. 
The  deforming  bodies  axe  schematically  illustrated  in  Figures  3(a)  and  (b),  which  also 
depict  the  initial  conditions  and  the  deformation  of  the  underlying  macromolecular 
network  structure.  The  various  loaded  and  relaxed  configurations  are  described  by 
their  appropriate  deformation  gradients  in  the  diagram.  For  the  case  of  the  initially 
isotropic  body,  the  initial  texture  of  the  undeformed  body,  Bo,  may  be  described  by  the 
identity  tensor,  I,  indicating  no  preferred  texture  in  the  body.  This  is  further  illustrated 
by  the  initially  isotropic  network  structure  shown  in  Figure  3(a).  In  other  words,  in  the 
model  of  Boyce,  et  al.,  the  initial  network  stretch  was  implicitly  assumed  to  be  identity. 
Earlier,  we  mentioned  that  the  initial  condition  on  the  athermal  shear  resistance  for  the 
material  in  the  annealed  state  is  given  by  «o  =  The  deformation  and  consequent 

kinematic  representation  of  the  initially  isotropic  body  was  discussed  earlier. 

For  the  corresponding  case  of  the  initially  oriented  body,  the  initial  texture  of  the 
body,  B0,  may  now  be  described  by  the  tensor  A’,  which  describes  the  network  preori¬ 
entation  of  the  body  in  terms  of  initial  stretch  ratios.  This  is  also  further  illustrated 
by  the  initially  oriented  or  stretched  network  structure  shown  in  Figure  3(b).  We 
wish  to  emphasize  here  that  this  initial  texture  was  not  necessarily  achieved  by  plastic 
working  of  the  material ,  but  could  have  been  otherwise  obtained,  c.g.,  by  a  spinning 
process.  The  body  is  now  deformed  via  F,  where  F  describes  the  applied  deformation 
gradient  relative  to  the  preoriented  state.  Returning  to  the  initially  textured  relaxed 
configuration  in  Figure  3b,  we  note  that  one  could  imagine  the  existence  of  a  "pre- 
reference”  configuration  possessing  am  isotropic  network,  such  that  both  the  state  and 
shape  of  the  textured  reference  configuration  Bo  could,  in  principle,  have  been  ob¬ 
tained  by  applying  the  (plastic  and  network)  stretch  A'  to  the  pre-reference  state  and 
configuration  The  deformation  of  a  material  element  relative  to  this  hypothetical  con- 
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figuration  may  be  described  by  FA'.  When  the  body  is  elastically  unloaded  without 
rotation,  the  relaxed  configuration,  relative  to  the  preoriented  state,  is  given  by  Fp. 
However,  the  deformation  of  the  relaxed  configuration  with  respect  to  our  hypothetical 
pre-reference  configuration  possessing  an  isotropic  network  structure,  would  be  given 
by  FPA'.  Therefore,  the  initial  texture  given  by  the  initial  network  stretch,  A*,  may  be 
operationally  interpretted  as  an  initial  plastic  stretch  tensor  (for  modelling  purposes). 
The  initial  condition  of  texture  is  incorporated  into  the  model  as  an  initial  plastic  de¬ 
formation  gradient.  Such  an  initial  condition  is  obtainable  in  a  glassy  polymer  by 
reduction  of  birefringence  data  (Kahar,  et  al.  [1978),  Ward  [1984],  Botto,  et  al.  [1984]). 
In  the  initially  isotropic  polymer,  the  back  stress  tensor  was  taken  to  be  coaxial  with 
the  left  plastic  stretch  tensor  because  the  latter  coincided  with  the  network  stretch.  In 
the  current  sense,  the  network  stretch  and  the  resulting  back  stress  is  a  result  of  textur¬ 
ing  which  occurred  both  prior  to  and  directly  from  the  currently  applied  deformation. 
Explicitly,  we  take  the  back  stress  tensor  to  be  coaxial  with  the  left  stretch  tensor  of 
the  product  FPA\  We  will  call  this  product  the  network  deformation  gradient,  F*. 
The  back  stress  would  be  obtained  as  follows: 

FN=FPA\  (9) 

Fn  =  V^R^,  polar  decomposition  (10) 

VN  =  R NTVNfLN,  eigen  value  decomposition  (11) 

Bi  =  BtiVf/X. i).  (12) 

Note  that  the  locking  network  stretch,  A*,,  is  that  of  the  isotropic  network  structure 
The  initial  value  for  the  athermal  shear  resistance  in  the  textured  polymer  must  be 
determined.  In  general,  the  method  of  preorientation  is  unknown.  Thus,  the  occur¬ 
rence  of  strain-softening  and/or  subsequent  physical  aging  is  also  unknown.  Therefore, 
currently,  the  initial  value  for  s  is  obtainable  by  conducting  a  stress-strain  test  and 
observing  if  strain-softening  occurs. 
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In  the  incorporation  of  preorientation  into  the  equations  of  the  constitutive  model 
given  earlier,  there  are  three  major  items  to  recall.  First,  the  deformation  gradient 
F  is  with  respect  to  the  initial  oriented  state.  Second,  the  back  stress  tensor  results 
from  all  contributions  to  molecular  alignment,  past  and  present.  In  other  words,  the 
back  stress  tensor  is  measured  with  respect  to  the  pre-reference  isotropic  state  of  the 
network.  Since  this  is  a  major  contribution  to  the  plastic  flow  resistance,  it  is  very 
important  to  get  a  good  measurement  of  the  initial  orientation,  A*.  The  third  item 
is  the  initial  value  of  the  athermal  shear  resistance,  which  may  or  may  not  be  that  of 
the  annealed  state.  Therefore,  the  constitutive  model  is  then  numerically  integrated 
and  incorporated  into  the  ABAQUS  [1984]  finite  element  code  as  done  in  Boyce,  tt 
al.  [1987b],  where  an  initial  condition  on  the  plastic  deformation  gradient  of  A*  is  now 
included. 

3  Effects  of  Preorientation  on  Material  Behavior 

Below,  we  will  quantitatively  examine  the  effects  of  texture  on  the  plastic  flow  of 
a  polymer.  During  homogeneous  deformation,  the  stress  at  which  active  plastic  flow 
begins  is  affected  by  both  the  degree  and  the  orientation  of  the  initial  texture.  Simulated 
tensile  tests  will  be  used  to  investigate  this  response.  Polymers  are  also  known  to  exhibit 
"shear  banding” ,  highly  localized  regions  of  straining,  during  plastic  flow.  Experiments 
( e.g .,  Brown  and  Ward  [1968],  Duckett  [1975],  Rider  and  Hargreaves  [1969])  have  shown 
that  the  direction  of  shear  banding  in  polymers  is  also  affected  by  both  the  amount  of 
texture  present  in  the  material  and  its  orientation  with  respect  to  the  principal  direction 
of  loading.  This  behavior  is  also  examined  with  our  constitutive  model  and  simulated 
tensile  tests.  The  numerical  results  are  compared  with  behavior  trends  reported  in  the 
literature. 
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3.1  “Homogeneous”  Behavior 


The  effect  of  orientation  on  the  initial  flow  behavior  of  a  glassy  polymer  is  now  ex¬ 
amined.  We  first  consider  some  experimental  results  from  tensile  tests  on  oriented 
polymers.  Secondly,  we  discuss  an  anisotropic  yield  criterion  which  has  been  applied 
to  such  polymers  in  the  past.  Finally,  using  the  constitutive  model  given  earlier,  we 
analyze  plastic  flow  in  an  oriented  glassy  polymer. 

3.1.1  Previous  Work 

Tensile  tests  on  thin  specimens  of  oriented  polymers  have  been  conducted  by  Rider 
and  Hargreaves  [1969,  1970]  on  amorphous  polyvinyl  chloride  (PVC),  by  Brown,  tt  al. 
[1968a,  b]  on  semi-crystalline  (nearly  amorphous)  polyethylene  terephthalate  (PET), 
and  by  Duckett,  tt  al.  [1972]  on  semi-crystalline  polypropylene  (PP).  These  investi¬ 
gators  have  taken  the  initial  texture  of  these  test  specimens  to  be  described  by  the 
principal  stretches:  Aj  =  A,  Aj  =  As  =  l/\/A,  where  A  is  called  the  initial  draw  rt  tio 
(IDR).  The  orientation  of  the  texture  of  a  specimen  when  tested  in  simple  tension  .s 
given  by  its  initial  draw  direction  (IDD)  which  is  the  angle,  9,  that  the  first  princip. 
stretch  makes  with  the  tensile  axis  (Figure  4).  Results  of  the  yield  stress  1  of  the  above 
materials  with  given  IDRs  as  a  function  of  IDD  when  tested  in  tension  (at  a  fixed 
displacement  rate)  are  shown  in  Figure  5.  These  materials  exhibit  a  significant  drop 
in  the  peak  stress  as  the  IDD  deviates  from  the  tensile  axis,  followed  by  a  gradual 
levelling  off  of  this  stress  at  an  IDD  of  approximately  9  *»  60°.  Brown,  tt  al.  and  Rider 
and  Hargreaves  then  determined,  from  their  experimental  data,  the  coefficients  for  an 
anisotropic  yield  criterion. 

A  criterion  for  yield  in  materials  which  have  been  preferentially  oriented,  t.g.  by 

deformation  processing,  such  that  they  retain  three  mutually  orthogonal  p lames  of  sym- 

1  Here,  we  ut  the  term  yield  itreee  to  describe  either  i)  the  peak  (trees  before  softening,  if  softening 
occurs;  or,  ii)  the  stress  at  the  point  where  the  tangents  to  the  elastic  and  hardening  portions  of 
stress-strain  curve  intersect. 
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metry,  has  been  proposed  by  Hill  [1950].  The  yield  surface,  which  reduces  to  the  Mises 
yield  surface  when  anisotropy  vanishes,  is  given  by: 


F{Tn  ~  TiS)*  +  G(Tis  -  Tlx)2  +  H{Tn  -  T„)2  +  2LT&  +  2 MTl  +  2NT*,  =  1;(13) 

where  Ti,  represent  the  stress  components;  the  subscripts  1,  2,  and  3  represent  the 
three  mutually  orthogonal  planes  of  symmetry,  where  the  direction  1  corresponds  to 
the  IDD  (Figure  4);  and  the  constants  F,  G,  H,  L,  M,  and  N  represent  the  square  of  the 
inverse  of  the  yield  strength  under  the  appropriate  normal  and  shear  loading  conditions. 
The  Hill  criterion  requires  identical  yield  points  in  tension  and  compression,  which  is 
generally  not  the  case  for  materials  exhibiting  a  Bauschinger  effect  and/or  pressure 
dependent  yield  —  including  polymers.  This  criterion  was  applied  to  the  PET  and 
PVC  experiments  discussed  above.  The  symmetry  of  the  pre-orienting  process,  where 
the  principal  stretch  ratios  were  as  given  earlier:  \i  =  A,  Aj  =  As  =  l/VX,  results 
in  the  equivalence  of  the  Hill  coefficients  G  and  H  for  this  case.  Due  to  the  loading 
in  these  experiments,  values  for  L  and  M  are  not  required.  Constanta  for  the  three 
remaining  Hill  coefficients  F,  G,  and  AT  were  obtained  from  the  data  of  the  0°,  45°,  and 
90°  IDD  tests.  The  resulting  curves  were  found  to  predict  the  tensile  yield  at  all  IDDs 
very  well.  These  investigators  noted  that  oriented  PET  and  PVC  behave  differently  in 
tension  and  compression,  exhibiting  a  Bauschinger  effect  as  well  as  pressure  dependent 
yield.  These  effects  are  not  accounted  for  in  the  Hill  Criterion.  Therefore,  both  Rider 
and  Hargreaves,  and  Brown,  et  al.  modified  the  criterion  to  include  a  scalar  internal 
stress,  or  back  stress,  6,  which  lumps  the  Bauschinger  and  pressure  effects  together. 
The  modified  criterion  then  reads  as  follows: 

jy(rn-6-r„),+/,(r„-rS3),+G(r»-r11+6),-h2Arr1,J+2Lr/s+2Mr*3  =  1.(14) 

Brown,  et  al.  also  conducted  tests  in  simple  shear  on  their  material.  The  modified 
criterion  was  found  to  give  good  agreement  with  the  shear  data  as  well  as  the  tensile 
data.  We  note  that  the  material  constants  F,  G,  N  and  5  are  for  a  given  initial  draw 
ratio,  as  well  as  at  a  given  strain  rate  and  temperature.  Once  these  material  constants 
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have  been  evaluated  for  the  desired  IDR,  yield  can  be  predicted  for  a  multiaxial  state 
of  stress  for  a  material  with  that  initial  texture.  If  one  is  relying  on  this  criterion  to 
predict  the  yield  behavior  as  a  function  of  the  initial  texture  in  the  material,  a  matrix 
of  the  material  parameters  F,G,H,L,M,Nt  and  b  must  be  assembled  by  conducting 
numerous  experiments  at  many  IDRs.  Furthermore,  these  coefficients  are  valid  only 
for  cases  of  pre-stretches  in  ratio  We  also  note  that  this  criterion  is  only  valid 

for  the  onset  of  yield  because  it  does  not  include  effects  of  continuing  deformation  on 
the  yield  parameters,  i.e.,  it  does  not  furnish  a  formalism  for  the  evolution  of  the  Hill 
coefficients  with  deformation.  Below,  we  discuss  the  prediction  of  tensile  yield  and  flow 
of  oriented  glassy  polymers  with  the  constitutive  model  of  this  paper.  An  important 
simplifying  result  is  that  the  model  is  based  on  material  properties  obtained  from  a  set 
of  experiments  on  the  material  in  the  isotropic  state  alone. 


3.1.2  Constitutive  Model  Results 

Here,  we  examine  the  plane  strain  tensile  yield  of  an  oriented  glassy  polymer*  where  the 
initial  texture  was  obtained,  for  example,  via  a  plane  strain  extrusion  process,  such  that 
the  texture  may  be  described  by  the  principal  stretch  ratios:  At  =  A,  Aj  =  1/A,  and 
A$  =  1,  where  A  is  the  IDR.  We  choose  to  consider  plane  strain  conditions  in  the  analyses 
in  order  to  maintain  a  two-dimensional  problem.  This  will  be  important  in  the  later 
simulation  of  shear  localization,  where,  had  we  chosen  to  model  the  tensile  testing  of 
thin  specimens,  the  problem  would  become  three-dimensional  once  localization  occurs. 
The  effects  that  the  amount  of  texture  and  the  orientation  of  the  texture  have  on 
the  initial  flow  response  of  the  material  are  inferred  by  considering  material  samples 
containing  various  principal  stretch  ratios  and  simulating  their  loading  at  a  number  of 
different  IDDs.  The  samples  of  material  are  assumed  to  have  been  given  enough  time 


3  For  illustrative  purposes,  the  material  constants  obtained  for  the  isotropic  PMMA  material  of  Boyce, 
et  aL  [1987a]  are  nsed  in  the  numerical  simulations  of  this  paper.  Analyses  are  conducted  at  room 
temperature,  where  we  recognise  that,  in  reality,  unoriented  PMMA  generally  does  not  shear  yield  at 
this  temperature,  but  erases,  and  fractures  prematurely. 
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to  have  aged  after  their  Initial  processing  and  prior  to  their  tensile  testing.  This  means 
that  any  softening  that  occurred  during  the  processing  which  produced  the  textured 
material  has  been  negated  (see  Boyce,  et  al.  [1987a]).  This  is  reflected  in  the  model 
by  using  an  appropriate  initial  condition  on  the  athermal  shear  resistance.  Here,  we 
assume  the  material  to  be  fully  aged,  so  the  value  for  s<>  of  the  initially  isotropic  material 
is  used.  Therefore,  the  material  again  softens.  The  orientation  and  texture  effects  are 
examined  by  comparing  the  stress  state  at  which  the  material  begins  to  flow  at  a  fixed 
temperature  and  strain  rate. 

For  the  case  of  plane  strain  tension,  we  compare  the  peak  flow  stress  prior  to 
softening  for  the  IDRs  of  A  =  1.0,  1.5,  2.0,  2.25,  and  2.5,  where  the  locking  stretch 
of  the  material  is  Ax,  =  3.0,  over  an  orientation  range  of  8  =  0"  to  8  =  90°  and  at  a 
constant  tensile  strain  rate  of  0.01  see-1  and  a  constant  temperature  of  25°C.  Before 
continuing,  we  note  here  that  homogeneous  deformation  does  not  actually  occur  in 
these  materials  after  yielding  because  of  the  large  tendency  for  deformation  to  localize. 
This  is  discussed  in  the  next  section.  We  conduct  this  analysis  in  order  to  obtain  the 
values  of  stress  at  which  plastic  flow  begins,  at  a  given  strain  rate  and  temperature, 
as  a  function  of  IDR  and  IDD.  The  results  are  shown  in  Figure  6.  The  impact  of 
orientation  is  found  to  increase  greatly  as  we  approach  initial  textures  near  the  locking 
regime,  i.e.  at  stretch  ratios  greater  than  2.0.  This,  of  course,  is  due  to  the  dramatic 
increase  in  the  back  stress  as  the  locking  stretch  is  approached.  We  also  note  that  these 
computed  curves  pass  through  the  isotropic  “yield”  point  at  an  orientation  8  less  than 
45s.  By  examining  the  flow  condition  for  plane  strain  tension  at  a  fixed  strain  rate  and 
temperature  and  neglecting  pressure  effects  (a  =  0),  we  see  that  the  back  stress  acts  to 
decrease  the  flow  stress  for  values  of  8  less  than  45s.  With  the  further  approximation 
F*  =  1,  we  obtain  r  to  be  [jCT'  -  B)  •  (T  -  B)]i,  where: 

'  0  0  o  ■ 

T  =  0  o  0  ;  (15) 

.0  0  f 


'  eos$  sin$  0  ]  [  61  0  0  1  f  coa ^  —sin6  0  ' 

B  =  —sin9  cos$  0  0  0  sin$  cosd  0  ;  (16) 

.0  0  lj[o  0  63  J  L  0  0  1. 

where  axe  the  principal  back  stress  components,  and  the  directions  1,  2,  3  and  the 
angle  9  were  shown  in  Figure  4.  The  flow  condition  is: 


n  0i  'Vo,# 
r  =  s  1  -  —  In—  «, 

As  ^ 


which  is  constant  for  a  given  strain  rate  and  temperature.  The  athermal  shear  strength, 
s,  is  taken  to  be  s0.  Incorporating  equations  (15)  and  (16)  into  (8)  gives: 


^<xJ  +  o(bx  -  6j)(sm*0  —  cos*9)  +  61s  +  6j*  +  =  2 r*.  (18) 

h 

For  the  case  of  0  =  45°,  the  tensile  flow  stress  is  a  =  2 yjr*  —  j(5iJ  +  6**  +  6j2).  There¬ 
fore,  at  am  IDD  of  45s,  the  flow  stress  will  be  less  than  the  isotropic  stress,  which  is 
given  by  a  =  2r.  This  decrease  in  flow  stress  is  accentuated  with  IDR  as  shown  in 
Figure  6.  The  results  for  the  larger  stretch  ratios  compare  well  with  the  trends  found 
in  the  experiments  discussed  earlier  and  shown  in  Figure  5. 

The  constitutive  model  can  be  used  to  examine  plastic  flow  given  any  general  multi- 
axial  loading  condition.  Yield  loci  cam  easily  be  obtained  for  any  initial  orientation  at 
any  fixed  strain  rate  and  temperature.  The  material  model  includes  the  effects  of  strain 
rate,  temperature,  pressure,  smd,  of  course,  the  additional  hardening  that  occurs  with 
continuing  plastic  deformation.  Therefore,  boundary  value  problems  containing  very 
general  and  inhomogeneous  loading  conditions  may  be  analyzed.  This  is  done  below 
by  investigating  the  effect  of  orientation  on  shear  banding  in  glassy  polymers. 


3.2  Localized  Flow 

Many  different  types  of  materials  are  known  to  experience  shear  localization  during 
plastic  flow.  The  occurrence  of  this  phenomenon  in  ductile  metals,  single  crystals, 
and  geological  materials  has  been  reviewed  by  Rice  [1976].  Shear  localization  is  a 
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phenomenon  which  is  of  great  interest  because  it  affects  such  a  broad  spectrum  of 
materials  and  can  have  a  marked  impact  on  all  subsequent  deformation  in  the  material, 
often  leading  to  fracture.  Solutions  determining  the  conditions  necessary  for  the  onset 
of  shear  localization  in  various  isotropic  material  models  have  been  obtained  for  the 
rate-independent  case  (Rudnicki  and  Rice  [1975],  Rice  [1976])  and  the  rate-dependent 
case  (Anand,  et  al.  [1987]).  Finite  element  analysis  may  also  be  used  to  predict  the 
direction  of  shear  banding  after  artificially  initiating  localization  with  either  a  material 
or  geometric  imperfection  (Tvergaard,  et  al.  [1981],  Needleman  and  Tvergaard  [1984]). 

3.2.1  Localization  in  Pre-Oriented  Polymers 

Experiments  on  pre-oriented  polymers  have  shown  that  shear  localization  is  affected 
by  both  the  degree  and  orientation  of  the  initial  texture  in  the  material.  This  effect 
has  been  examined  experimentally  by  Duckett  [1975]  and  Brown  and  Ward  [1968]  in 
the  fully  oriented  semi-crystalline  (nearly  amorphous)  polymer  PET  below  its  0,;  by 
Duckett,  et  al.  [1972]  in  the  highly  oriented  semi-crystalline  polymer  PP  above  its 
Qt;  and  by  Rider  and  Hargreaves  [1969]  in  the  amorphous  polymer  PVC.  In  these 


experiments,  textured  polymers  were  tested  at  various  EDDs  (and  various  IDRs  for  the 
PVC  case)  in  simple  tension.  The  resulting  shear  band  direction  (SBD)  was  recorded. 
This  direction  is  given  by  the  angle  i//  shown  in  Figure  4.  The  experiments  on  the 
semi-crystalline  polymer  sheets  showed  the  SBD  to  be  closer  to  the  tensile  axis  than 
the  IDD  except  for  IDDs  of  8  <  16°,  in  which  cases  shear  banding  was  difficult  to 
observe,  t.e.  homogeneous  deformation  was  obtained.  Rider  and  Hargreaves’  plots  of 
the  SBD  as  a  function  of  IDD  of  PVC  at  various  draw  ratios  are  shown  in  Figure  7. 
The  data  is  from  tensile  tests  on  thin  sheets  which  have  been  oriented  such  that  the 
principal  stretch  ratios  are  given  by:  A;  =  A,  Aj  =  Aj  =  l/y/X.  This  experimental 
data  indicate  that  the  SBD  depends  on  the  specific  material,  the  initial  draw  ratio,  and 
the  orientation.  As  the  IDR  increases,  its  impact  on  the  direction  of  shear  localization 
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becomes  stronger,  especially  as  the  locking  stretch  ratio  is  neared.  This  shows  that 
the  SBD,  is  dependent  on  both  material  and  IDR,  i.e.,  the  effect  of  IDR  on  the 
material  response  is  relative  to  the  locking  stretch  of  the  specific  material.  We  can  see 
these  effects  in  the  amorphous  PVC  data  where  the  IDR  has  a  very  large  impact  on  the 
SBD  as  a  function  of  orientation,  0.  Figure  7  shows  that  at  a  low  IDR,  the  SBD  is  not 
significantly  impacted  by  the  orientation  $  of  the  IDD,  and  the  0  of  the  SBD  becomes 
a  strong  function  of  0  as  the  IDR  increases  and  approaches  the  effective  locking  stretch 
of  the  material. 

These  investigators  (Duckett  [1975],  Brown  and  Ward  [1968],  Duckett,  tt  al.  [1972]) 
have  been  able  to  analytically  predict  the  direction  of  shear  banding  in  these  polymers 
with  reasonable  accuracy  by  using  Hill’s  anisotropic  theory  of  yield  (equation  5)  and  a 
normality  (associated)  flow  rule.  The  resulting  analysis  predicts  two  directions  for  the 
formation  of  a  shear  band.  These  predicted  directions  agreed  well  with  the  experimental 
data  and  were  further  improved  upon  by  using  the  Hill  criterion  as  modified  with  a 
scalar  back  stress,  which  was  described  earlier.  We  emphasize  that  using  this  method 
requires  experimentally  obtaining  the  Hill  coefficients  and  the  scalar  back  stress  at  each 
draw  ratio. 

Once  the  material  has  localized  and  a  prominent  shear  band  forms,  the  material 
within  the  band  reorients  toward  the  tensile  axis.  This  phenomenon  has  been  exper¬ 
imentally  observed  in  PET  (Brown  and  Ward  [1968]),  PP  (Duckett,  tt  al.  [1972]), 
and  PVC  (Rider  and  Hargreaves  [1970])  by  examining  the  deformed  material  with  a 
polarizing  microscope.  The  extinction  direction  corresponds  to  the  maximum  refrac¬ 
tive  index  which  is  parallel  to  the  (local)  principal  stretch  direction.  The  extinction 
direction  in  the  band  was  found  to  be  different  from  that  of  the  bulk  material.  The 
extinction  direction  within  the  band  rotated  from  the  IDD  towards  the  tensile  axis  by 
an  angle  (0-9),  where  9  is  the  angle  that  the  principal  stretch  ratio  within  the  band 
currently  makes  with  the  tensile  axis.  This  data  indicated  that  the  material  within 
the  band  reorients  towards  the  loading  axis.  Brown  and  W?*d  had  additional  results 
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showing  that  the  greater  the  amount  of  applied  tensile  deformation,  the  greater  the 
amount  of  reorientation.  This  was  more  clearly  demonstrated  by  Rider  and  Hargreaves 
in  their  results  on  PVC,  which  are  shown  in  the  plots  of  Figure  8.  These  plots  show 
the  orientation,  $,  of  the  material  within  the  band  with  respect  to  the  tensile  axis  as 
a  function  of  the  applied  nominal  stretch  ratio.  The  nominal  stretch  is  determined 
from  the  total  applied  displacement  and  the  initial  gage  length.  We  note  that  the 
local  stretch  in  the  shear  band  would  be  much  greater  due  to  the  localization  of  the 
deformation  in  the  region.  The  initial  orientation  at  unit  stretch  is  the  angle  of  the 
EDD.  These  results  indicate  that  the  material  within  the  band  continues  to  reorient 
as  it  is  stretched  until  the  orientation  direction  is  parallel  with  the  tensile  axis.  Rider 
and  Hargreaves’  results  also  demonstrate  that  the  reorientation  occurs  faster,  i.e.  at 
lower  nominal  stretch  ratioe,  when  the  IDR  is  smaller.  This,  of  course,  is  due  to  the 
fact  that  there  is  less  re-orienting  to  accomplish.  Our  analysis  below  will  examine  the 
phenomenon  of  reorientation  within  the  shear  band. 

3.2.2  Numerical  Simulation  of  Localization 

The  case  of  a  plane  strain  tensile  test  of  oriented  PMMA  at  room  temperature  was 
analyzed  here  using  the  constitutive  model  of  Boyce,  et  al.  [1987a]  with  appropriate 
initial  values  for  A*  and  «.  A  point  of  principal  interest  was  the  effect  of  material 
orientation  on  the  direction  of  shear  banding  in  glassy  polymers.  A  relatively  coarse 
mesh  of  4-node  plane  strain  elements  (ABAQUS  type  CPE4)  was  used  in  the  finite 
element  analysis.  This  element  is  the  constant  dilation  4-node  element  of  Nagtegaal, 
et  al.  [1974],  where  the  shear  and  pressure  terms  are  integrated  separately  to  prevent 
“locking”  of  the  mesh.  We  note  that  a  sharper  shear  band  could  be  picked  up  with 
a  finer  mesh  and/or  a  mesh  of  quadrilateral  elements  composed  of  constant  strain 
triangle  (CST)  elements  (Needleman  and  Tvergaard  [1984]).  The  element  size  defines 
the  minimum  possible  thickness  of  the  band.  Therefore,  an  analysis  with  a  finer  mesh 
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of  4-node  elements  could  result  in  a  narrower  shear  band.  The  shear  band  will  form 
along  element  boundaries.  The  CST  mesh  formation  provides  a  greater  freedom  for 
the  deformation  to  localize  because  of  the  larger  number  and  direction  of  element 
boundaries.  Also,  because  these  are  constant  strain  elements,  large  strains  can  occur 
in  an  element  without  necessarily  affecting  its  neighboring  elements.  This,  of  course, 
is  precisely  the  physical  response  of  shear  banding.  If  one  is  interested  in  obtaining  a 
narrow  band,  and  the  direction  of  the  band  is  known  a  priori,  then  one  could  tailor 
a  CST  mesh  to  accurately  pick  up  the  localized  straining.  However,  the  purpose  of 
this  analysis  was  to  compute  the  direction  of  the  shear  banding.  Therefore,  the  coarse 
unbiased  mesh  of  quadrilateral  elements  shown  in  Figure  8  was  chosen,  although  it  may 
result  in  an  artificially  broad  shear  band. 

The  simulation  to  be  discussed  is  for  the  plane  strain  tensile  testing  of  a  pre-oriented 
glassy  polymer,  where  the  plane  of  straining  during  the  pre-orienting  process  and  the 
subsequent  testing  are  identical.  The  problem  becomes  asymmetric  with  respect  to 
the  tensile  axis  when  the  principal  direction  of  orientation  is  aligned  at  some  angle 
with  respect  to  the  tensile  axis.  The  asymmetry  of  the  problem  which  results  from 
material  pre-orientation  permits  a  model  of  one-half  of  the  test  specimen  to  be  used 
with  appropriate  boundary  conditions.  This  condition  is  enforced  by  constraining  the 
nodes  on  the  bottom  of  the  mesh  on  either  side  of  the  specimen  center  to  displace 
in  equal  and  opposite  x  and  y  directions.  Due  to  the  material  orientation,  shear  band 
formation  in  the  central  portion  of  the  specimen  may  result  in  rotation  of  the  specimen. 
This  was  permitted  to  occur  by  constraining  the  nodes  along  the  top  of  the  mesh  to 
remain  on  a  straight  line,  but  by  giving  them  the  freedom  to  contract  as  well  as  pivot 
about  the  center  node  of  this  surface.  This  is  equivalent  to  loading  the  specimen  in  a 
pinned  grip.  The  center  node  of  the  bottom  surface  of  the  mesh  is  constrained  in  the 
x-  and  y-directions.  The  center  node  of  the  top  surface  of  the  mesh  is  constrained  in 
the  x-direction  only  and  is  displaced  in  the  y-direction  to  the  desired  extension.  The 
half-length  of  the  specimen  was  chosen  to  be  four  times  the  width  in  order  to  allow  the 
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shear  band  to  grow  unhindered  along  a  straight  path  until  reaching  a  free  surface. 

The  shear  band  is  initiated  in  the  center  of  the  specimen  by  making  the  two  bottom 
center  elements  of  the  mesh  softer  than  the  remainder  of  the  mesh.  This  was  done 
by  assigning  these  elements  a  lower  value  for  the  material  property  of  the  athermal 
shear  resistance,  so  =  105 MPa,  than  the  rest  of  the  elements  of  the  mesh,  which  were 
assigned  so  =  112MPo.  This  caused  a  shear  band  to  emanate  from  this  point  once 
the  material  began  to  flow.  Different  perturbation  values  of  so  <  112 MPa  were  tried 
on  these  two  elements  in  a  test  case  to  ascertain  its  effect  on  the  resulting  shear  band 
direction.  It  was  found  to  have  no  effect. 

The  analysis  was  conducted  on  a  mathematical  model  with  the  same  material  con¬ 
stants  as  determined  by  Boyce,  et  al.  [1987a]  for  PMMA.  The  material  was  taken  to 
be  oriented  to  varying  levels  of  texture  such  that  the  principal  values  of  the  initial 
textures  were  given  by:  Ai  =  A,  Aj  =  1/A,  Aj  =  1,  where  A  is  the  initial  draw  ratio,  and 
Aa  coincides  with  the  plane  strain  direction  of  the  simulation.  The  material  was  then 
“tested”  in  plane  strain  tension  at  the  initial  draw  ratios  of  A  =  2.0,  2.25,  and  2.50, 
and  the  initial  draw  directions  of  $  =  0°,  30®,  45®,  60®,  and  90®.  The  direction  of  shear 
banding  was  obtained  by  examining  contours  of  the  plastic  strain  rate  where  regions 
of  intense  plastic  straining  indicate  the  development  of  a  shear  band.  The  angle  that 
these  regions  make  with  the  tensile  axis  gives  the  direction  of  banding.  Two  bands  of 
concentrated  straining  initially  form  in  each  sample.  These  bands  are  at  approximately 
90®  to  one  another  with  possibly  a  few  degrees  difference  (Anand  and  Spitzig  [1982]) 
due  to  the  pressure  dependent  yield  of  the  polymer  which  is  included  in  the  constitu¬ 
tive  model.  The  angle,  0,  at  which  the  shear  band  formed  with  respect  to  the  tensile 
axis  is  depicted  in  Figure  10,  where  it  is  shown  to  change  with  both  initial  draw  ratio 
and  direction.  These  numerical  results  demonstrate  that  as  the  IDR  approaches  the 
locking  regime,  where  the  locking  stretch  of  PMMA  is  A*  =  3.0,  the  dependence  of 
the  shear  band  direction  on  orientation  becomes  more  acute.  This  follows  the  trends 
found  in  experiments  discussed  earlier,  which  showed  that  the  more  highly  textured 
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the  polymer,  the  greater  is  the  dependence  of  shear  band  direction  on  orientation. 

Another  interesting  result  of  the  numerical  analysis  is  illustrated  in  Figure  11.  This 
figure  shows  the  contours  of  plastic  strain  rate  as  deformation  progresses  for  the  cases 
where  the  IDR  is  initially  oriented  with  the  tensile  axis.  For  the  high  draw  ratio  of 
2.5,  we  see  that  the  material  at  first  begins  to  localize  into  symmetric  45°  shear  bands; 
however,  the  straining  soon  becomes  uniform  throughout  the  test  specimen.  This  was 
also  observed  in  the  experiments  on  highly  oriented  semi-crystalline  PET  where  clearly 
defined  shear  bands  were  not  observed  at  orientations  of  9  <  15°.  Duckett,  et  al.  did 
not  comment  on  any  test  they  may  have  conducted  at  0°  orientation  on  polypropylene. 
Rider  and  Hargreaves  did  not  obtain  distinct  banding  for  A  =  2.0  and  9  =  0°  or  for 
A  =  3.3  and  9  —  0",  15®  in  their  PVC  samples.  It  appears  that  they  did  not  test  the 
highest  draw  ratio  at  these  orientations  for  this  reason. 

The  absence  of  shear  banding  in  these  experimental  cases  as  well  as  the  numerical 
simulations  at  0®  orientation  can  be  explained  by  analogy  with  a  Considgre  stability 
analysis.  In  the  isotropic  material,  necking  occurs  in  the  Considgre  construction  at 
the  point  when  the  slope  of  the  stress-strain  curve  equals  the  stress,  <x,  divided  by 
the  stretch,  A,  or  ^  =  |.  Stability  is  regained  and  cold-drawing  begins  at  a  second 
Considgre  point.  The  highly  preoriented  material  which  is  aligned  with  the  tensile  axis 
is  essentially  at  a  state  past  the  second  Considgre  point  and  in  the  middle  of  a  stable 
cold-drawing  mode,  which  even  the  artificial  initiation  of  localization  cannot  supersede. 
Referring  to  the  stress-stretch  curve  for  isotropic  PMMA,  our  material  of  IDR=2.5  was 
past  this  second  Considgre  point  whereas  the  lower  IDR  of  2.25  was  before  but  quite 
near  this  point.  Therefore,  we  can  expect  the  material  with  IDR=2.5  to  continue  in 
the  stable  drawing  mode  and  the  material  with  IDR=2.25  to  require  a  greater  amount 
of  applied  strain  before  localization  disappears  and  stable  cold-drawing  takes  over. 
This  is  exactly  what  was  predicted  with  the  model  as  illustrated  in  Figure  11,  which 
demonstrates  that  a  larger  applied  stretch  was  needed  before  homogeneous  drawing 
occurred  in  the  material  with  an  IDR  of  2.25  than  in  the  material  with  an  IDR  of  2.5. 
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Our  model  correctly  predicts  this  continued  drawing. 

It  is  also  of  interest  to  observe  the  development  of  a  shear  band  in  those  cases  where 
the  IDR  is  aligned  asymmetrically  about  the  tensile  axis.  The  case  of  am  IDR  of  A=2.5 
and  an  IDD  of  8=  30“  was  chosen  to  examine  this  development.  Figures  12  and  13 
depict  the  contour  plots  of  plastic  shear  strain  rate  and  plots  of  the  deformed  mesh 
as  deformation  progresses  for  this  case.  For  purposes  of  illustration,  the  asymmetry 
conditions  were  used  to  depict  both  halves  of  the  specimen.  The  plots  of  the  deforming 
mesh  show  the  deformation  to  be  concentrated  in  a  shear  band.  The  material  asymme¬ 
try  also  causes  a  slight  shearing  of  the  bulk  of  the  specimen.  The  plastic  shear  strain 
rate  contours  illustrate  the  increase  in  intensity,  with  deformation,  of  straining  within 
the  band,  where  the  plastic  shear  strain  rate  is  as  much  as  five  times  greater  than  the 
applied  nominal  tensile  strain  rate  of  O.Olsee-1.  The  rate  contours  for  the  case  of  0 
=  30°  also  show  the  broadening  of  the  band  where  the  region  of  highest  shear  strain 
rate  moves  away  from  the  center  of  the  specimen.  This  indicates  that  the  material 
near  the  specimen  center,  where  the  shear  band  initially  formed,  is  now  essentially 
“locked”.  Therefore,  active  plastic  deformation  spreads  up  (and  down)  the  specimen, 
thus,  broadening  the  band.  We  also  observe  two  concentrated  bands  of  intense  strain¬ 
ing  in  the  initial  stages  of  deformation.  However,  as  deformation  progresses,  one  band 
disappears,  and  we  eventually  form  a  single  distinct  shear  band  through  the  first  and 
third  quadrants  of  the  specimen,  which  also  contains  the  IDD.  This  is  in  agreement 
with  experimental  observations.  However,  the  second  band  which  appeared  in  the  ini¬ 
tial  stages  of  our  numerical  analysis  is,  on  occassion,  the  shear  band  which  forms  in 
the  experiments.  We  note  that  direct  co  mparisons  between  our  plane  strain  numerical 
results  and  the  experimental  results  on  thin  tensile  specimens  cannot  be  made  due  to 
the  different  constraints  of  these  tests,  where  thinning  of  the  bands  is  permitted  in  the 
experiments. 

We  were  also  able  to  predict  the  reorientation  of  the  material  within  the  band.  The 
amount  of  reorientation  was  found  to  depend  on  the  applied  deformation.  Figure  14 


I.VI.V 


.‘k,  ****4i‘ 


■Mil 


depicts  the  computed  reorientation  of  the  material  within  the  band  for  the  case  of 


an  IDR  of  2.5  and  an  IDD  of  30°  as  a  function  of  the  applied  stretch.  We  note  here 
that  the  applied  stretch  (nominal)  is  relative  to  the  total  specimen  length  of  8.0  units, 
whereas  the  active  plastic  deformation  occurs  over  a  gauge  length  of  approximately  2.6 
units.  We  see  the  material  reorients  toward  the  tensile  axis  by  a  greater  amount  as  the 
applied  tensile  stretch  increases.  These  results  are  in  accord  with  those  found  in  the 


experiments  discussed  earlier. 


4  Concluding  Remarks 


In  this  paper,  the  effect  of  preorientation  on  the  inelastic  deformation  behavior  of 
glassy  polymers  was  examined.  This  effect  was  incorporated  into  the  constitutive  model 
of  Boyce,  Parks,  and  Argon  [1987a]  via  the  use  of  appropriate  initial  conditions  on 
internal  state  variables.  This  method  of  modelling  texture  in  a  material  takes  full 
advantage  of  the  concept  of  internal  state  variables.  The  model  was  subsequently 
used  in  conjunction  with  the  finite  element  method  to  analyze  the  homogeneous  and 
localized  flow  of  oriented  glassy  polymers.  Comparing  the  results  of  our  numerical 
simulations  to  trends  found  in  similar  experiments,  we  found  that  our  model  correctly 
assesses  the  initial  flow  behavior,  as  well  as  several  complicated  responses  resulting  from 
the  inhomogeneous  process  of  shear  localization.  These  include  the  direction  of  shear 
banding  as  a  function  of  initial  orientation  conditions,  material  reorientation  within 
the  band,  and  the  conditions  which  result  in  continued  cold  drawing  of  the  specimen. 
The  numerical  analyses  did  not  require  any  further  experimental  material  property 
identification  beyond  the  isotropic  material  properties  and  the  initial  texture. 

We  have  shown  that  the  model  in  its  current  form  is  capable  of  successfully  capturing 
many  of  the  major  effects  of  initial  texture  on  inelastic  deformation  in  glassy  polymers. 
However,  in  order  to  obtain  a  complete  constitutive  description,  many  items  still  need 
to  be  addressed.  One  such  item,  important  for  both  initially  isotropic  and  oriented 


polymers,  is  the  effect  of  the  physical  aging  process  on  plastic  flow.  Temperature  and 
strain  rate  have  also  been  observed  to  affect  the  strain  hardening  response  of  glassy 
polymers  (Hope,  et  al.  [1980],  Botto,  et  al.  [1987]).  This  effect  is  considered  to  be 
the  result  of  a  relaxation  of  the  network  structure  (Raha  and  Bowden  [1972],  Kahar, 
et  al.  [1978],  Boyce  [1986],  Botto,  et  al.  [1987]).  Additionally,  extreme  terminal 
hardening  has  been  observed  in  crystallizable  polymers  such  as  PET,  whose  initial 
structure  was  amorphous  (Ward  [1984]].  This  degree  of  hardening  cannot  be  entirely 
explained  by  the  backstress  of  the  current  formulation  and  could  be  of  major  importance 
in  the  response  of  highly  oriented  polymers.  These  features,  which  require  further 
experimental  quantification  and  modelling,  are  currently  under  investigation. 
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0,  INITIAL  DRAW  DIRECTION  (DEGREES) 

jure  5.  Nominal  tensile  yield  stress  as  a  function  of 
initial  draw  direction,  for  highly  oriented  PET  and  PP, 
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6,  INITIAL  DRAW  DIRECTION  (DEGREES) 

gure  6.  Computed  peak  flow  stress  as  a  function  of  0,  initial  draw 
rection,  for  various  initial  draw  ratios  (IDR)  for  PMMA  at  ©=25°C. 
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Figure  9.  Finite  element  model  of  one-half  of  a  plane  strain  tensile  test  specimen. 
Boundary  conditions  consist  of  constraining  nodes  on  either  side  of  the  specimen 
center  along  the  bottom  surface  of  the  mesh  to  displace  in  equal  and  opposite  x- 
and  y-directions.  The  nodes  along  the  top  surface  are  constrained  to  remain  on 
a  straight  line.  The  top  center  node  is  constrained  in  the  x-direction. 
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within  the  shear  band  as  a  function  of  nominal  stretch  for 


an  IDR  of  2.5  and  IDD  of  30°  for  PMMA. 
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